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Abstract. We report on a 3600 s spectrum of GRB 021004 obtained with the Nordic Optical Telecope on La Palma 
10.71 hours after the burst. We identify absorption lines from five systems at redshifts 1.3806, 1.6039, 2.2983, 
2.3230, and 2.3292. In addition we find an emission line which, if due to Lya from the host galaxy, gives a redshift 
of 2.3351. The nearest absorber is blueshifted by 530 km s~^ with respect to this line, consistent with shifts seen 
in Damped Lya and Lyman-Break galaxies at similar redshifts. The emission line flux is 2.46 ± 0.50 x 10^^^ 
erg s~^ cm~^. Some of the absorption systems are "line-locked", an effect often seen in QSO absorption systems 
believed to originate close to the QSO central engine. 
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1. Introduction 

More than 30 Optical Afterglows (OAs) to Gamma-Ray 
Bursts (GRBs) have been detected to date. Absorption 
lines due to metal enriched gas in the host galaxy of the 
GRB have been detected in all cases where a spectrum 
with a good signal-to-noise ratio of the OA has been ob- 
tained. Prior to GRB 021004 the redshift has been de- 
termined from such detected metal absorption for about 
a dozen OAs. Absorption line spectra of GRBs have pri- 
marily been used to constrain the redshift of the burst but 
also to gain insight into the nature of the burst, the gas 
phase of its host galaxy, and of its large scale environment 
(e.g. Vreeswijk et al. 2001; Mirabal et al. 2002). 

Furthermore, it is interesting to consider how the GRB 
absorption line systems compare to the well-studied QSO 
absorption line systems, and to examine whether the phys- 
ical conditions (e.g. column densities, metallicities, ioniza- 
tions states, and kinematics) in the absorbers detected in 
GRB afterglow spectra are similar to any of the various 
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classes of QSO absorbers (e.g. Jensen et al. 2001; Savagho 
et al. 2002a; Salamanca et al. 2002). 

In this Letter we present optical spectroscopy of the af- 
terglow of GRB 021004. GRB 021004 was detected by the 
HETE-II satellite on October 4.5043 2002 UT (Shirasaki 
et al. 2002). An optical afterglow was reported already 
10 minutes after the burst (Fox 2002). The first optical 
spectroscopy reported two absorption systems at z=1.38 
and 1.60 (Fox et al. 2002), which were later confirmed 
by Eracleous et al. (2002), Anupama et al. (2002) and 
Castander et al. (2002). Lya absorption at z=2.323 was 
first reported by Chornock & Filippenko (2002). The 
presence of this system was confirmed by Salamanca et 
al. (2002b), Savaglio et al. (2002b), Sahu et al. (2002) 
and Castro-Tirado et al. (2002). Its lightcurve showed 
strong deviations from the usual power-law decay indi- 
cating structure in the surrounding medium (Lazzati et 
al. 2002; Holland et al. in prep). 

2. Observations and data reduction 

2.1. Observations 

The spectrum of GRB 021004 presented in Fig. 1 
was obtained with the 2.56-m Nordic Optical Telescope 
(NOT) using the Andalucia Faint Object Spectrograph 
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Fig. 1. Normalised spectrum of GRB 021004. Absorption lines are marked and numbered as in Table 1. 



(ALFOSC). The detector was a 2048^ pixels thinned Loral 
CCD with a pixel scale of 0'/189. We used a grism cover- 
ing 3850A- 6850 A and a slit of 1.3 arcscc providing a 
resolution of 6.5 A and a sampling of 1.5 A per pixel. 

The spectroscopic observations consisted of a single 
3600 s exposure started at Oct. 4.9505 (10.71 hrs af- 
ter the burst). Photometric observations obtained on the 
same and the following nights at the NOT and elsewhere 
made it possible to closely follow the GRB light-curve. 
The brightness of the OA at the start of the exposure was 
R=18.48±0.03 and B=19.15±0.10. During the one-hour 
exposure the OA faded by 0.08 magnitudes. These broad 
band magnitudes were used for the flux calibration. 

2.2. Data reduction and absorption line-list 

Standard techniques were used for bias, dark and flat 
field corrections. For the optimal ID extraction, and for 
the line-search, measurement, identification, and fitting 
we used a code originally developed for 2D spectral PSF 
fitting (M0ller 2000) and for QSO spectral analysis (for 
details see M0ller & Kjaergaard 1992). 

In Table 1 we list centroid (vacuum corrected) and 
observed equivalent width of all absorption lines found 
above our 5a detection limit. Several of the listed lines 
are complex blends, and for those we provide the total ob- 
served quantities with no attempt at this point to deblend 
them. Inferred redshifts are listed only for the single un- 
blended lines. The absorption lines are marked above the 
normalised spectrum shown in Fig. 1. 

We have identified five absorption systems at redshifts 
1.3806, 1.6039, 2.2983, 2.3230, and 2.3292 and we shall in 
what follows refer to them as systems A, B, C, D, and E. 
The identification of those systems was carried out inde- 
pendently of the early GCNGs, but is seen to agree well 
with the circulars summarised in Sect. 1. 

For systems C, D, and E our spectrum covers Lya. No 
line redwards of the Lya line of system E is unidentified, 
and we conclude that line No 4 in Table 1 marks the onset 
of the Lyman forest. System E is therefore most likely as- 
sociated with the GRB itself, and could mark the redshift 
of the GRB host galaxy. There are however a few caveats 
and we shall return to this question in Sect. 4. 



Table 1. Absorption lines in GRB 021004 



No. 


Avac 


Wobs 


aw 


line ID (system) 


^ahs 




(A) 


(A) 


(A) 






1 


3981.75 


7.10 


0.63 


complex (Fig. 2) 




2 


4011.90 


13.11 


0.54 


complex (Fig. 2) 




3 


4026.32 


2.46 


0.31 






4 


4041.43 


14.18 


0.47 


complex (Fig. 2) 




5 


4352.06 


2.79 


0.30 


Alii (B) 


1.6048 


6 


4398.21 


2.95 


0.35 


Cii (C) 


2.2957 


7 


4438.65 


2.57 


0.34 


C II (D-l-E) 




8 


4596.82 


1.26 


0.21 


Si IV (C) 


2.2982 


9 


4634.89 


7.97 


0.29 


Si IV (C+D+E) 




10 


4667.64 


4.95 


0.28 


Si IV (D+E) 




11 


5108.88 


5.96 


0.26 


Civ (C) 




12 


5152.16 


19.56 


0.26 


C IV (D-l-E) 




13 


5511.21 


0.84 


0.20 


Alii (C) 


2.2986 


14 


5562.56 


1.99 


0.19 


Alii (E) 


2.3293 


"15 


5582.76 


2.41 


? 


Fell (A) 




16 


5652.16 


1.32 


0.19 


Fell (A) 


1.3804 


17 


5672.26 


2.38 


0.17 


Fell (A) 


1.3805 


18 


6102.42 


1.64 


0.23 


Fell (B) 


1.6032 


19 


6158.50 


2.41 


0.23 


Fell (A) 


1.3809 


20 


6190.76 


7.58 


0.33 


complex (Fig. 2) 




21 


6658.57 


4.19 


0.22 


Mgii (A) 


1.3812 


22 


6675.02 


3.83 


0.21 


Mgii (A) 


1.3809 


'23 


6734.36 


0.84 


? 


Fell (B) 




24 


6768.53 


2.00 


0.18 


Fell (B) 


1.6031 



" Line 15 is influenced by stron, 
Line 23 is strongly influenced 



g sky subtraction residuals, 
by a cosmic. 



Most of the absorption lines are blended, making it im- 
possible to obtain a good redshift directly from line cen- 
troiding. The redshifts given above were therefore derived 
through line profile fitting. Detailed fits to the blends are 
shown in Fig. 2, and the redshifts required to fit each 
species are listed in Table 2. 

The lines are heavily saturated and the resolution low 
enough that the line profiles are dominated by the reso- 
lution. Therefore the column densities used in the fits are 
mostly too poorly constrained to be useful. Nevertheless, 
from the Lya line profiles we can set a strict upper limit 
of 1.1 X 10^° cm~^ on the Hi column density. The total 
Hi column actually used in the fit of systems C, D, and 
E was 5.7 X lO^^ cm'^. 
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Table 2. Absorption redshift constraints from line fitting 





Fell 


Civ 


All! 


Si IV 


Total 


A 


1.3806 








1.3806 ± 0.0005 


B 


1.6030 




1.6048 




1.6039 ± 0.0009 


C 




2.2980 


2.2990 


2.2982 


2.2983 ± 0.0004 


D 




2.3220 




2.3240 


2.3230 ± 0.0010 


E 




2.3290 


2.3290 


2.3296 


2.3292 ± 0.0004 



s ^ cm ^ identical to other reported measurements to 
within la. The flux is comparable to the Lya flux of the 
GRB 000926 host galaxy (Fynbo et al. 2002). The inferred 
redshift (2.3351) is 530 km s~^ higher than that of the 
E system. At the time of writing the limit on the host 
galaxy magnitude is B>24, providing a lower limit on the 
observed emission line equivalent width of 170 A. 
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Fig. 2. Profile fits and line identifications. 



3. Lya in emission 

Inspection of the top panel of Fig. 2 shows that what is 
mostly a very good fit of the Lya absorption lines of sys- 
tems C, D, and E, becomes very poor on the red shoulder 
of the last line. Subtraction of the fit (Fig. 3) reveals that 
the poor fit is caused by the presence of an emission line 
shortly redwards of the Lya-E line. As suggested in several 
GCN Circulars (Chornock & Filippcnko 2002; Salamanca 
et al. 2002b; Djorgovski et al. 2002; Castro-Tirado et al. 
2002), this line is likely Lya emission from the host galaxy. 
The line (Fig. 3, right panel) is unresolved, has its cen- 
troid at 4054.4 A, and a fiux of 2.46 ± 0.50 x lO-^^ erg 





Fig. 3. Left: Section of the spectrum with Lya absorp- 
tion lines. Note the poor fit to the extreme right wing of 
the second line. Right: Residual after subtraction of the 
fit shown to the left. A narrow Lya emission line is now 
clearly visible. The vertical scale in both plots is in units 
of 10"^^ erg s"^ cm-2 A'^ 



4. Discussion 

4.1. Properties of Lya absorption and emission lines 

A shift of a few hundred km s~^ between redshifts de- 
termined from Lya absorption and emission is commonly 
seen for DLA galaxies (M0ller et al. 2002) and for Lyman- 
Break galaxies (e.g. Adelberger at al. 2002, their Figures 
2 and 3). Two effects arc thought to contribute to this ob- 
served velocity shift: if the absorbing gas is in a galactic 
wind moving towards us, that will blueshift the absorber 
and at the same time cause the resonantly scattered Lya 
photons to become redshifted in order to escape the cloud. 
A spread in absorption/emission redshifts as large as that 
seen in systems C-E and host Lya emission (3300 km s~^) 
has never been seen in Lyman-Break Galaxy spectra, and 
in fact equation 2 from Adelberger et al. (2002) shows 
that the velocity shift of the Lya emission line should 
be less than 200 km s~^ for an emission line rest equiv- 
alent width larger than 50 A. Therefore the large veloc- 
ity range spanned here cannot be explained in terms of 
normal host galaxy properties, and we conclude that the 
absorbers must be in the local environment of the GRB. 

4.2. Tlie nature of the GRB absorption systems 

Any absorption seen in the spectra of high-rcdshift QSOs 
may be classified as belonging to one of three basic cate- 
gories (Weymann et al. 1979). The "intervening systems" 
arc cosmologically distributed and can have any redshift 
lower than that of the source, "local cluster systems" will 
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have the same redshift as the source to within ± the veloc- 
ity dispersion of the local cluster or filament, and "ejected 
systems" (physically closest) can be radiatively acceler- 
ated to large velocities. 

In QSOs the ejected systems are often seen as BALs 
(P-Cygni type broad absorption line systems) and/or as 
line- locked systems (Foltz et al. 1987). Ejection velocities 
as high as 0.1c have been reported (Vilkoviskij & Irwin 
2001), and the systems arc found to be highly ionized and 
to have extremely high metallicities (Savaglio et al. 1994; 
M0ller et al. 1994). The high ionization is easily under- 
stood because of the intense UV flux of the QSO. The 
line-locking is explained as an effect of cloudlets being 
accelerated via absorption at a given transition until its 
wavelength falls into the shadow of another line in a cloud 
in front of it (Vilkoviskij et al. 1999). 

In this scheme we would classify systems A and B as 
intervening systems, while systems C, D and E display 
a surprising similarity to the ejected systems of QSOs. 
They all have strong absorption of highly ionized CIV 
and Si IV but no detectable Sill absorption, they have 
high column densities compared to stellar winds, and at 
intermediate resolution they show evidence for line-locking 
in C IV, Si III/Lya and possibly Si IV (see also Savaglio et 
al. 2002b; Salamanca et al. 2002b). For confirmation of the 
line-locking higher resolution spectra (already obtained) 
are needed. 

It is difficult to understand why there is this resem- 
blance between systems C-E and ejected systems of QSOs. 
We concluded above that the systems arise in the local 
environment of the GRB. They may therefore be related 
either directly or indirectly to the GRB or the GRB pro- 
genitor itself, or they may be caused by an unrelated phe- 
nomenon which just happens to be sharing the same vol- 
ume of space. Firstly we would consider it unlikely that 
the systems are directly related to the burst itself. So soon 
after the burst any material ejected simultaneous with 
the GRB itself should undergo significant changes on rel- 
atively short timesealcs, yet there is no evidence for such 
changes between our spectrum and that of Matheson et 
al. 2002. Secondly one may ask if stellar winds could cre- 
ate such signatures. Radiatively driven winds with termi- 
nal velocities of up to several 1000 km s~^ are ubiquitous 
for very hot stars such as O stars and WR stars (Lucy 
& Abbott 1994; Kudritzki 2002) but these winds cannot 
reach the required column densities. Also, presumably for 
the same reason, actual line-locking has never been ob- 
served in stellar winds. 

We cannot exclude that the GRB is located close to 
the inner region of a recently deceased QSO, especially 
at Z'^2, where the QSO density is high. However, given 
the existing evidence that GRBs seem to be related to 
the deaths of massive stars it is more likely that the GRB 
progenitor was a massive, hot star and that this star radia- 
tively drove the fast outfiow now observed in absorption 
against the light of the afterglow. Other massive stars and 
supernova explosions could be contributing to the wind if 
the progenitor was located in a compact star-forming re- 



gion similar to that of GRB 980425 (Fynbo et al. 2001). 
A similar scenario was suggested for GRB 971214 by Ahn 
(2000). 

Alternatively, as suggested by Lazzati et al. (2002), in 
the supra-nova scenario (Vietri & Stella 1998) the wind 
and the clumpy surrounding medium could be the result 
of a supernova predating the GRB by several years. 
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